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We discuss several points that may help to clarify some questions that remain about the anomaly 
puzzle in supersymmetric theories. In particular, we consider a general Af = 1 supersymmetric 
Yang-Mills theory. The anomaly puzzle concerns the question of whether there is a consistent way 
to put the 7?-current and the stress tensor in a single supercurrent, even though in the classical 
theory they are in the same supermultiplet. As is well known, the classically conserved supercurrent 
bifurcates into two supercurrents having different anomalies in the quantum regime. The most 
interesting result we obtain is an explicit expression for the lowest component of one of the two 
supercurrents in 4-dimensional spacetime, namely the supercurrent that has the energy-momentum 
tensor as one of its components. This expression for the lowest component is an energy- dependent 
linear combination of two chiral currents, which itself does not correspond to a classically conserved 
£NJ . chiral current. The lowest component of the other supercurrent, namely, the 7?-current, satisfies 

the Adler-Bardeen theorem. The lowest component of the first supercurrent has an anomaly that 
we show is consistent with the anomaly of the trace of the energy-momentum tensor. Therefore, 
we conclude that there is no consistent way to put the 7?-current and the stress tensor in a single 
supercurrent in the quantized theory. We also discuss and try to clarify some technical points in the 
derivations of the two-supercurrents in the literature. These latter points concern the significance 
of infrared contributions to the NSVZ /3-function and the role of the equations of motion in deriving 
the two supercurrents. 

PACS numbers: ll.10.Hi, 11.30.Pb, 11.30.Rd 

I. INTRODUCTION 

m 
> 

' The anomaly puzzle in Af — 1 supersymmetric gauge theories is well known. A real superfield, J^, called the 
\ supercurrent can be constructed [l[ and is classically conserved. The lowest component of this superfield is the R- 
C*") ■ current. The other components of are related to the supersymmetry current 3 a{i (where a is a two-component 
spinor index that labels the generators of the supersymmetry) and the stress tensor -d^ through linear transformations. 
. This construction is related to the fact that these symmetries are elements of the superconformal algebra. 

The anomaly puzzle arises as follows. In an M = 1 SYM (supersymmetric Yang-Mills) theory, the R-symmetry, 
which is just a chiral U(l) symmetry (denoted later as U(1)r) has an anomaly. This chiral anomaly is proportional 
to the topological invariant, F^F^, and can be expressed in an operator equation. One can try to generalize the 
operator equation of this anomaly of the i?-symmetry to a supersymmetric form involving [2|-|7j]. However, this 
attempt led to an apparent contradiction. On the one hand, the anomaly of R-symmetry is known to be exactly 
of one- loop order because of the Adler-Bardeen theorem @, Q. On the other hand, the trace of the stress tensor, 
which is another component of D a \J a & should be proportional to the /3- function (because the trace is a measure of 
the breaking of scale invariance). These two components of F) a J a a should be proportional to the same factor, which 
would seem to imply that the /3-function is exactly of one loo p o rder. However, explicit perturbative calculations 
show that there are higher order corrections to the /3-function [10| . Note that there are some subtleties about this 
formulation of the anomaly puzzle, which we shall discuss in more detail later. 

In Grisaru et al, [TP . Il2j , a solution to the anomaly puzzle is given by showing that there are actually two different 
supercurrents J^. Let us call those two different supercurrents in 4-dimensional spacetime, and J^ 2 \. They are 

the same classically (meaning at tree level). One of them, J^ 1 ' has the R-current as its lowest component, but the 
higher components are no longer the supersymmetry current and stress tensor. The anomalous non-conservation of 
this supercurrent is proportional to the one- loop /3-function. The other supercurrent, J^ 2 \, has the supersymmetry 
current and stress tensor as its components and has an anomaly proportional the exact /3-function. In Ensign et al 
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13 1 , th ey consider J\f = 1 supersymmetric gauge theories including matter fields and extend the construction done in 
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12j of the two supercurrents to the case that includes matter. 
In Shifman and Vainshtein [3] , they argued for a different solution to the anomaly puzzle. They considered the co- 
efficient in front of the W 2 term in the Wilsonian effective action and showed that its running, obtained by integrating 
out higher momentum modes, is only of one-loop order. (This result is in agreement with the nonrenormalization the- 
orem, as formulated for example in [la].) On the other hand, the physical coupling constant, defined from the physical 
amplitudes one measures in experiments, includes higher-order contributions. Motivated by this, they proposed an 
operator anomaly equation with a one-loop coefficient. They showed that the /3-function of the physical coupling is 
obtained when one takes the matrix element of their anomaly equation. They assumed that their single operator 
supercurrent contains as its lowest component the R-current and as another component the stress tensor. However, 
we show below that this is not the case, by proving that the supercurrent, having the energy- momentum tensor 

as one of its components does not, in fact, have the R-current as its lowest component. To avoid any ambiguity, we 
mention that we are using the term i?-current (and i?-symmetry) to describe the U{1) current (denoted by R^) that 
transforms the gaugino A, the matter scalar A and the matter spinor tp according to the charge ratios of 1 : | : — |. 
We find that J^ 1 ^ has as its lowest component a current which is an energy- dependent linear combination of the 
R-current and the Konishi current. This linear combination, which we refer to as i?^, is not a chiral current. This 
favors the idea of two supercurrents proposed in [H1 - fl3l [. The explicit expression for the lowest component of this 
supercurrent had not been written earlier to our knowledge. 

The anomaly equation for SYM with matter fields, as given in [3], has a term 7_D 2 ($e v $) (where $ is a chiral 
superfield) that is responsible for the anomalous dimensions of the matter fields. This term does not appear in [l3| 
because they assume that external fields are on-shell. As we shall see, it is the existence of this term that implies that 
the lowest component of R' is not the i?-current but a mixing (with coupling constant dependent coefficients) of the 
.R-current and the Konishi current. We perform an explicit calculation, which is not in the literature, to obtain the 
mixing. We perform the calculation using both component fields and the supersymmetric background field method. 
The results we obtain from either method agree and give the r yD 2 ($e v <f>) term. 

The Adler-Bardeen theorem implies that the chiral current i? M has an anomaly of one-loop order, but does not 
imply that R' should have a one-loop anomaly. We also find that the difference between R' and 7? M is manifest in 
a very clear way at the infrared fixed point, where R'^ becomes a non-anomalous symmetry current that is a linear 
combination of R^ and the Konishi current. We use the supersymmetric QCD model as an example to show this 
behavior of the operator R' . 

Although the two-supercurrent scenario appears to be the correct solution to the anomaly puzzle, there are some 
technical issues in their construction that we discuss and attempt to clarify. In [l2( , the equations of motion (EoM) 
are applied with the assumption that they vanish (up to contact terms). However, if one uses the expectation values 
of the various operators, as given in then the EoM would seem to have nonvanishing expectation values. We 
show that this apparent inconsistency is resolved when one takes into account the non-local contributions. After 
that, the expectation values of the bare operators are consistent with the application of the EoM. In particular, 
the expectation value, (V a j7ad), of the unrenormalized operator V a J a a vanishes as required to by the EoM. More 
explicitly, the non-local contribution to (V q J7qq) is opposite in sign to the local contribution, which is proportional to 
an e dimensional operator , and the two contributions add up to zero in the limit that e — > 0, i.e., in 4 dimensions. As 
a result, (V q 1 7 Q q) does vanish. Then, when we use the renormalization procedure of [TT| - fT3| . in which the contribution 
proportional to an e dimensional operator is removed by renormalization, the non-local contribution indeed gives the 
correct one-loop anomaly. This correct one-loop anomaly was obtained in [TT| - [l~3l j . They did not explicitly discuss the 
role played by the non-local contributions in their derivation, so the discussion of those terms here may help clarify 
the consistency of the construction of the two supercurrents. 

Finally, we comment on the question of whether the higher-order terms in the /3-function are the result of contribu- 



tions coming from infrared modes of the fields. In [14[, they show that the higher-order terms in the /3-function come 



from the infrared modes. A different way of obtaining the same /3-function is given in [l6| . In the latter method, the 
coupling constant receives its higher-order corrections from the Jacobian appearing when one rescales the measure 
[l7l [Hj], and as they mention in the method does not appear to depend on the infrared modes. By changing 
the UV cutoff in the Wilson effective action, we show that the momentum modes above any arbitrary finite non-zero 
scale do not give a significant contribution to the Jacobian from which the multi-loop corrections to the /3-function 
are obtained. This shows that the method used by [16| does indeed depend on the infrared modes. 

In section II, we review some basic ideas about the supercurrent and the anomaly puzzle. The supercurrent is 
discussed in more detail in the appendix. In section III, possible solutions to the anomaly puzzle in the literature 



1 The calculation is performed using dimensional reduction and the dimension is 4 — e with e > 0. 
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are reviewed and remaining problems are discussed. In section IV, we perform an explicit calculation to show that 
the operator R' in the same supermultiplet as the supersymmetry current has exactly the properties of what the 
anomaly equation in [T3] predicts but it generates a U(l) transformation different from the i?-symmetry. As a 
result, this supcrfield should be identified as M and not as M (in the notation defined above). First we do 
the calculation using component fields. Then in section IVA, we obtain the same result using the supersymmetric 
background field method. In section IVB, we analyze the properties of the current R' at the non-trivial infrared fixed 
point of supersymmetric QCD. We show that R' does have the charge ratios to be a non-anomalous current and thus 
corresponds to a true symmetry at the fixed point, as it should. In section V, we discuss the role of non-local terms 
in obtaining the expectation value of the equation of motion and show how such terms enter into the construction of 
the two supercurrents. In section VI, we show that the calculations of the /3-function done by Shifman and Vainshtein 
[l4| and by Arkani-Hamed and Murayama [l6[ both depend on the infrared modes. 



II. REVIEW OF ANOMALY PUZZLE 



First of all, let us review the basics of the supercurrent and the anomaly puzzle. Af = 1 gauge theory is described 
by a Lagrangian, 

£ = » s / d 2 9TvW 2 + H.c. 
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T(R) 

= / d 2 9TrW 2 +H.c. (1) 

V J 

T(R) denotes one half of the Dynkin index for the representation R, Ti(T a T b ) = T(R)8 ab . The superfield W a = W£T a 
in components is 2 

W a = -\D 2 e- v D a e v = -i{X a + iO a D + 6? f a p + i99V a6i \ & ), 
o 

where f a p is the field strength in the spinor coordinate, f a p = ~\{(J tl cr' / ) a pF lll , = —a^F^. The vector superfield V 
in the Wess-Zumino gauge is, V = -29 a 6i v a a + 2i999\- 2i999\ + 8066(D + id^v^). The integration over Grassmann 
numbers is defined by J 9 2 d 2 9 — 2. Recall that the anomaly puzzle can be stated in terms of the absence of a 
supersymmetric anomaly equation. Such a possible equation is described by a supercurrent J^, which is a superfield 
and can be defined as [14| 

Jaa = ~^Tr[e v W a e~ v W«] = -—^—^[e v W a e~ v W & ]. (2) 
9 9 T{R) 

Generally, the components of the supercurrent superfield are related to the i?-current R^, the supercurrent J af2 and 
the stress tensor "d^ respectively, 

Jad = C a& + {e^X0c t c t + H.C.} + 2e^^T a& ^-^{eJ $ id^C 1 a l + H.C.} 

+{\^M a& + H.c} + {^6 2 6% a& + H.c} + U 2 9 2 D a6t . (3) 

C M , Xn an d T ^v are related to i? M , J M and as we shall see in the appendix. 
For the supercurrent defined by Eq.([2|), we have, 

4 

C aa = Rad = oTr(A a Ad), (4) 

g 2 

The 99 component of J a a, (O corresponds to the stress tensor and the exterior derivative of R^. Note that 
is not really the stress tensor as we shall see in the appendix. However, the trace r M M is proportional to that of the 
stress tensor . So the operator also gives the trace anomaly. 



2 We mostly use the conventions in Wess and Bagger 19] including the choice of cr^ matrices and superderivatives. Our conventions differ 
only with regard to the normalization of the vector superfield and the integration of Grassmann variables. 
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As usual, this R-symmetry is broken at the quantum level because it is a chiral U(l) symmetry. The anomaly 
equation is, 

^R„ = ~^F^F afa/ , (5) 

as follows from the Adler-Bardeen theorem Q. T{G) is the T{R) of the adjoint representation. 
One can lift Eq. (|5j) to the supersymmetric form, 

d^J* 1 = l -CTr[D 2 W 2 - D 2 W\ (6) 

where C is some coefficient to be determined. The lowest component of ([6]) is the chiral anomaly equation, Eq.©. 
Equivalently, we have, 

D & J a& = CD a TrW 2 . (7) 

The real part of the 9 component D a J a a in Eq.© corresponds to t m m while the imaginary part corresponds 
to — <9 M i? M . This matches the component of D a TrW 2 , whose real and imaginary parts are —e a pTr (FF) and 
— e Q( gTr(FF) respectively. 

To be consistent with the prediction of the Adler-Bardeen theorem, both sides should be bare operator and C 
should be of one-loop just like Eq.©. However, to get the correct trace anomaly, /3-function, which has higher-loop 
contributions should appear on the right hand side. So C has to be proportional to /3. Now we get the anomaly puzzle. 
At least, this is how this puzzle was originally stated. There are quite a few subtleties as we shall see. 

The situation becomes more complicated when matter is introduced. The Lagrangian becomes, 

£=-L / d 2 6TvW 2 +H.c. +i f d i 6j2¥e v <i> f . (8) 
9 J J f 

where $ f are chiral supcrfields describing matter. The supercurrent is defined as 

Jad = -^Tr[e v W a e~ v W & ] + ^ ^ ¥ (v A e y V Q - e v D«V a + vj5 a e v ) (9) 
9 * f 

where covariant derivative is introduced V a $/ = e~ v D a e v $> /. The i?-current has the form, 

= 4Tr(Ac7„A) - i {^f a ^S ~ 2iA f ^ M A}) , (10) 
9 ^ f V 

where Af is the scalar component of the chiral superficld and tpf is the spinor component. 

With the introduction of matter, there is another U(l) symmetry <&/ — > e la Qf. The corresponding current is the 
so-called Konishi current (denoted by Rj[). This symmetry is certainly chiral and its anomaly, the Konishi anomaly 
is given by, 

D 2 jf = D 2 (<t>fe v <t> f ) = ^ILttW 2 . (11) 
Equivalently, we can define a superfield Q a & as in Eq. pl?]) , which has R? as its lowest component. 



III. POSSIBLE SOLUTIONS TO THE PUZZLE 

Alternately, we can use Wilson effective action to describe the anomaly puzzle [lij]. In this scenario, the theory has 
a large but finite cutoff. The Wilson effective action at scale A is denoted by SV(A). Higher momentum modes can 
be integrated out to provide the running of the coupling constant. It can be shown that the new SV(A — <5A) obtained 
by this renormalization group flow will only have a one-loop correction to the coefficient of TrVK 2 in the Lagrangian 
©. This agrees with the conclusion based on the non-renormalization theorem [l5|. However, this result appears to 
be in contradiction to the multi-loop (3 function. Note that the coefficient of the Wilson effective action (at scale A) 
can be related to the 1PI amplitude with an infrared cutoff A (see e.g. [20]). As noted by Shifman and Vainshtein, 
the absence of the infrared modes is the reason for the absence of multi-loop corrections in the /3-function. Shifman 
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and Vainshtein distinguish between the physical coupling constant and the corresponding coefficient in the Wilson 
effective action Syy- The latter is renormalized only at one- loop level as predicted by the nonrenormalization theorem. 
On the other hand, the physical coupling can be obtained by evaluating the matrix elements (or the effective action) . 
To do that, all the infrared modes have to be included and the higher-order corrections emerge. 

Shifman and Vainshtein then proceed to propose that a single supercurrent can contain both the stress tensor 
and the i?-current. The anomaly equation for this bare supercurrent J a a is of the form of Eq.Q with a one- loop 
coefficient C (see e.g. Eq.(19) in [14j). 



2 

n a t ■ — - n 

— U aa — 3 



2 „ 
= -D a 
3 



16tt 2 

/ 



(12) 



The {go) is the one-loop /3-function and (3g 1 \go)/go 1S a go-independent number. Note that the operators in the 
equation are bare operators. To obtain the physical coupling constant one needs to take the matrix elements of the 
operators on the right hand side. The matrix of element of W 2 is shown to have finite multi-loop contribution that 
exactly reproduce the correct /3-function. More explicitly, when the operators on the right are expressed in terms of 
renormalized operators, Eq. (|12l) becomes the anomaly equation that has the correct multi-loop /3-function, 



9 

where [ ] indicates renormalized operators, and 



D«J a * = 2 Da {^l[TrW 2 } - |5>We V */)]}, (13) 

/ 



g3 3T(G)-Z f T{Rf){l-if) 
^ = '16^ I-tW/2, • (14) 

The term proportional to 7^ comes from the second term in Eq. (|12l) because of the Konishi anomaly (jlip . More 
explicitly, the contribution (proportional to 7/) from the second term in Eq. (|12j) to the /3-function follows when the 
operators on the right hand side are diagonalized. 

Note that the trace of the stress tensor, should be equal to ^ j3 a {M)O a {M), in which O a {M) are renormalized 
operators (at scale M). Moreover, [2l| the coefficients of these operators O a {M) can be considered as /3- functions (with 
g a {M) as variables) only when the operators O a {M) are "orthonormal" at scale M. More explicitly, the operators 
O a {M) are chosen so that the corresponding matrix element 3 of every coupling constant g a only receives contribution 
from a single operator (on the right hand side of the trace anomaly) and the matrix element should be exactly unity 
(up to some power of M's). Only in this case we can take the coefficients of the operators on the right hand side of 
the trace anomaly to be the ^-functions. 

However, there are some subtleties about Ea. (TT2l that imply a contradiction with the Adler-Bardeen theorem. Let 
us look at it more carefully. For the example of pure SYM, the 9 2 component of an operator W 2 in fact has an 
imaginary part equal to (where ao = g 2 /47r), 



~FF - 0„(A<7"A) = ~FF - 4 7rao ^ff l : 



where we used the fact that the second term on the left hand side is proportional to <9 M i? M (see Eq.Q). After this term 
is moved to the left side in Eg. (TT2"|) . it is clear that Eq. p^| does not reproduce the Adler-Bardeen theorem; namely 
that the anomaly of the R-current is no longer proportional (with a coupling-constant-independent proportionality 
factor) to the topological term FF as in the non-supersymmetry gauge theory. In other words, Eq.© no longer holds 
as an operator equation of bare operators. Moreover, unlike it was previously claimed in the literature [22j, Eq. (|13[) 
does not agree with the Adler-Bardeen theorem. Even if Eg- fTS"]) is not obtained from the one- loop equation ([12")) 
but is taken as the starting point, what appears on the right are just renormalized operators and can not be moved 
to the left side, which only contains bare operators. On the other hand, it has been shown [23j that if the lowest 
component of the supercurrent on the left of Eg. (1121) is taken as a renormalized operator, which is different from the 



3 Without the insertion, this matrix element corresponds to a certain amplitude that defines the coupling constant at scale M. 
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bare i?-current by a multiplicative renormalization factor, correct anomaly equations (for both trace anomaly and 
chiral anomaly) can be obtained. 

There is another way to show the inconsistency between Eq.(fT2|) and the Adler-Bardeen theorem. Together with 



;ncy 

the proposed expectation value of W 2 ((46) in [lj]) 



2tt 



t2 

' ext > 



equation (TT^j) predicts a nonvanishing expectation value for the bare chiral current (d^R^) at two-loop level. More 
explicitly, the two-loop value is ^ a times the one-loop value. This conclusion however, is in contradiction with the 
combination of the Adler-Bardeen theorem and the proposed expectation value of FF ((57) in [Hj]), 



(FF) = (FF). 



T(G)a 



which implies a two-loop expectation value being T ^ G ' a times the one-loop value. Note that the Adler-Bardeen 
theorem states that Eq. ([5]) is an operator equation of bare operators and the expectation values of both sides should 
have the same quantum corrections 4 . Such an agreement is confirmed up to two- loop in M. 

So eventually we have no choice but to construct two supercurrents. One of them, jWj, has the i?-current as its 
lowest component, but does not have the stress tensor among its components, while the other J 1 ^ ^ has the stress 
tensor but not the i?-current. As a result, there is no reason to have a single operator equation to describe both 
chiral anomaly and trace anomaly. The construction of two supercurrents using the background field method and 
dimensional reduction is first proposed by Grisaru et al [ll| f° r the pure SYM and is further developed by Ensign 
and Mahanthappa [13] for the coupled SYM. As we shall see later, there is an inconsistency in their calculation. 
However, we show by careful calculation that their results for the two currents are indeed correct. 

Let us briefly review their results. In this approach, two different renormalized currents (both superfields) are 
defined. Each satisfies an anomaly equation with the anomaly expressed in terms of renormalized operators. One of 
the anomaly equations is similar to Eq. (|13j) with the renormalized coupling constant and operators on the right hand 
side. The other has an one-loop coefficient for the W 2 term, which agrees with the Adler-Bardeen theorem 5 . The 
explicit form is (Eq.(3.16) in [13j after a change in convention), 

V Qd [J Qd ] = [[V a W' 3 V a W ] - [V & W^ & W^) . (15) 

Both currents are renormalized operators whose expectation values are finite. 

We shall see how these two supercurrents are constructed. For later convenience, we will give the operators that 
are involved, 

W a6l = -^Tr[e- v W & e v W a ] (16) 

K a& = W a& - ^Tr piVaW" - *T Q £VdW^ (17) 
P a& = i($e v V Qti $ - $V a(i e y $) (18) 
Q a * = ~\V a ,V&K*e v $). (19) 

To be consistent with the expressions given above in section II, we use the covariant derivatives in the gauge chiral 

representation. Note that T a a introduced in [l2j is the e-dimensional projection of the gauge connection and is gauge 
covariant under the K gauge transformation (and invariant under the A gauge transformation). Now in the gauge 
chiral representation, it is covariant under the A gauge transformation. The expectation values of various operators 
are given by (3.3) (3.4) (3.5) in [13]. They are obtained by the background field method and dimensional reduction. 



4 Hi gher order quantum corrections to FF in QED are discussed in [23 |. 

5 In and |13|| . only the divergent contribution to the expectation value is considered and therefore the Adler-Bardeen theorem implies 
the absence of any two-loop contribution to the coefficient C. 
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The dimension is d — 4 — 2e. Some of the expectation values given below are different from from those in [13|. This 
is due to different conventions. For example, one of the equation given by (3.4) in [13J is, 

WJ (2) =0xC + -, (20) 

The superscript e denotes (renormalized) external fields. The superscript r denotes renormalized fields ($ = Z^ 2 $ r ) 
with the field strength Z$ given by, 

Z $ = 1 + 2 9 —C(R) + r 4 ( 1 1 ^ ' ^ 1 ^' T>\ r r* i m , 

e 

and Zy given by 



Z* = l + 2*-C(R) +9 {-£--) (ST(G)C(R) + C(R)T(R) + 2C(R) 2 ), (21) 



o 2 o 4 
= 1 + y[3T(G) - T(R)] + [3T(G) 2 - T(G)T(R) - 2C(R)T(R)]. 

However, {Q 1 '^) 1 ^ a re not the two- loop expectation values of the Q r aa . Instead, it is the two-loop expectation value 
of the operator renormalized to one-loop order. In this paper, we use the symbol (O)^ for the expectation value 
of an operator O without subtracting any subdivergence due to renormalization of this operator. In this convention, 
Eq. (|20| is expressed as 

(Q r adl + $1 (2C(R)Q r a& + T(R)K r a .))W = x K* a & + . . . , (22) 

where C(R) is the quadratic Casimir operator of representation R. Note that only the two-loop contribution propor- 
tional to K a a is evaluated and the rest is unknown. With field strength renormalization Z$ (Q a a = Zq>Q r aAl ) given 
by Eq.jUl) and {K r a£t ) {l) = 0, Eq.([22]) can be rewritten as, 

(Q aa ) {2) =0xKZ & + .... 

Let us work with SQED, in which the corrections to the expectation value of FF start at two- loop and the corrections 
to the right hand side of the Adler-Bardeen theorem start at three-loop. Naively, one can speculate that the U(l) 
current in Q a & satisfies the Adler-Bardeen theorem in the sense that there is no anomaly at two-loop level. Note that 
(Qaa)^ 1 ' is nonvanishing and leads to a nonvanishing expectation value of <9 M C^ (C$ being the lowest component of 

Q ^ . . i-i 

However, Q aci is not the correct superfield containing the anomalous U(l) current. In the approach used by [13j, 

has to be renormalized and the anomaly is described by a renormalized operator [Qa]. However, in the usual anomaly 
calculation of non-supersymmetric gauge theories (with matter), the expectation value of a bare chiral current <9 M j^ is 

proportional to FF. Anyway, if we ignore this difference and just apply the equations of motion on those bare fields 
from which [Q^] is constructed, the correct anomaly equation follows. 

Following Eq.(3.8) in fl3j |. we can find out the relationship between the bare operators and the renormalized 
operators as, 

W aa = [W aa }~^(T(R)[W aa ]^T(G)[K]-3C(R)[P aa ]-C(R)[Q aa ]) 



£(3T(G)T(R) + C(R)T(R)) + 9 1{ 1 -T{G)T{R) + \ 



P aa = [P aa ] ~ — (3C(R)[P aa ] + C{R)[Q aa ] - T{R)[W a& ]) 
e 

9 -{ZT{G)T(R) + C(R)T(R)) - ^(T(G)T(R) - C(R)T(R))] [K a6l ] 

2 

Qaa = [Q a a\~ 9 —T{R)[K a6 ] (23) 

Renormalized operators like [W Q q] are defined to have the expectation values of the background fields. With Eq. 
we get to the conclusion that a current 

\J act — act Pact ~t~ TiQaai 
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has no anomaly at two-loop level because of the lack of g 4 [K aa } in (J„„). However, this bare operator J aa contains 
neither nor 9^ . Instead, two different renormalized operators need to be constructed to describe the two anomalies 
(trace and chiral). One of them is the supercurrent (J'*- 2 -' M in our notation), 

Jaa — ^^aa Paa ~l~ aa 0(e), (^4) 

where the 0(e) terms do not affect the /3-function. The other is the Adler-Bardeen current (J^ 1 -^ in our notation), 

[J a6l ] = [W aa ] + [Pad] + g[Qod]- 

As shown in [l3j], with the use of the equations of motion, the desired anomaly equations (fTBf and (|T5|) can be obtained. 

A major problem this approach has is about the use of equation of motion. For example, the trace anomaly is 
described by Eq. ([M|) . The left hand side is a renormalized operator in the sense that expectation value of the operator 
[Jaa] is given by the background fields. There is no way this expectation value can give what appears on the right 
hand side under a derivative V". In fact, the right hand side is obtained by taking the expectation value of another 
operator W 2 which is obtained by the equation of motion of K aa ■ Kaa is in the definition of the renormalized operator 
[Jaa] (see Eq. (f!M)) ). The rest of Jaa, the o per ator W aa + Paa + \Qaa gives no anomalous contribution because of the 
EoM. This is very confusing. As shown in jl3j ]. the two sides of an equation of motion generally do not have the same 
expectation values. For example, V a W aa = (for pure SYM) follows from the equations of motion but apparently 
V a (yV a a) 7^ following from Eq.([23f (see also Eq.(3.3) in [13(). This problem will be further discussed in section V 
and a possible solution will be proposed. 



IV. SUPERPARTNER OF THE TRACE ANOMALY 



Before we move on to talk about the solution to the problem about the equation of motion in the construction of 
the two supercurrents, let us give a supporting argument for this approach. As explained in section III, the one-loop 
anomaly equation Eq. (TT2l implies that the operator R'^, defined as the lowest component of does not satisfy 
the Adler-Bardeen theorem. So it is unlikely that R'^ is the the i?-current R^ 6 . Here we try to use some explicit 
calculation to show that for a general coupling g, this operator R' is a mixing of the current i? M and the Konishi 
current R^ {R^ being the U(l) current in Q a a)- This result clearly supports the approach of two supercurrents. It 
also gives a clear physical interpretation of the lowest component of J aa (the one containing i?^„), which is not given 
before. This validity of this interpretation is particularly clear at the infrared fixed point where the superconformal 
symmetry is restored. At this point, the charges of those fields A, A and if) under the i?-symmetry are different from 
their classical values. This new U{1) symmetry is also a classical symmetry whose current R'^ is a linear combination 
of R^ and the Konishi current Rj[. The latter assigns charge +1 to both A and ip. Moreover, as we shall explain, this 
property of R'^ agrees with the last term in the anomaly equation Eq. (I12[) , which is actually not obtained in [13j . So 
the J aa in the Shifman-Vainshtein scenario should be identified as the supercurrent with z9 M „. 

To study R'^, we compute the Green's functions of this operator (or rather d^R 1 ^) and various other fields. More 
explicitly, we compute the Green's functions with an insertion of the operator obtained via supersymmetry transfor- 
mation of the gamma trace (see below) of the supersymmetry current. This operator has a term d^R'^ according to 
the superconformal algebra, 

{S a , Qp} = 4M a p - 2iDe a0 - 3R'e a p. (25) 

where M a p are the Lorentz generators and S a is the generator corresponding to the gamma trace of J^. We will now 
show that the contact terms of the Green's function, which are the changes of the other fields under the transformation 
generated by R'^, can be described by the transformation of a combination of the original P-symmetry and the Konishi 
U(l) symmetry. 

Let us start with the computation of cr M J M . The gamma trace of the supersymmetry current in the Wess-Zumino 
model is, 

a% = -ia^{Xf, + W v Xv) = 3^x M = -2\f2a»d^A -> -2s/2a»V^A (26) 



Note that we define ij-current as the (7(1) current associated with the (anomalous) symmetry that transforms the gaugino A, the matter 
scalar A and the matter spinor ip according to the charge ratios of 1 : j : — which is determined by the classical supercurrent. 
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In the last step, we include the effect of the gauge field by covariantizing the derivative. Eq. ([26|) corresponds to the 
contribution from the matter multiplet to the gamma trace of the full supersymmetry current. 

However, Eq. ([26|) does not vanish on-shell because of the interaction with the gauge field. Let us consider SQED 
for simplicity. The equation of motion of ip is, 



So the gamma trace in 4d is 



fflu^ip = -V2ieA*X. 



a »J = -2V2a^VaipA - 4eAA*X, 



or equivalently, we should add — ^ea^XAA* to the definition of m Eq. ([26|) . Of course, this term can also be 
obtained from explicit calculation (from a Xa^D term in the supercurrent of the gauge multiplet). In 4 + e dimension, 
the gamma trace becomes, 



V v Aa w a^ + ^ v d v (Aip) - ^[--ea^XAA* + a^a" (--ea^XAA*)] 



V2{(2 + e)V„Aat l i>+~ 



(4 + e) 



(2 + e)] d^(Aa^)} - (4 + e)(l + ^)eAA*X 



^eV^Aa^ - ^eAa^V^j - (e + j)eAA*X 



V2 



eD tl Aa tt ip-eeAA*X. 



(27) 



The supersymmetry transformation (parametrized by £ Q ) of the gamma trace is (note that the indices in a v £ are 



-(■D tl ipa'"ip - 2V2e AipX + 2iV fl A*V»A) 
3 



2V2 



eA*i>X 



eV2eAi/jX + e^-A*tpX + 



3A/2e _ , , 
e— — A*ipX + E.c. 



(28) 



We use the equation of motion in the middle step. This is justified because 

(V^ipa^ + V2eAX)if), 

is the the counting operator [25j, which has a finite expectation value. In the last step, we just keep the real part 
that is needed (and drop the imaginary term V il A*V^A). Note that this is because the imaginary part should be 
proportional to the trace anomaly following the superconformal algebra. 

The Green's functions with an insertion of 6^(<J^ J M ) can be evaluated. Alternately, the same results can be obtained 
by calculating the expectation value of S^^a^J^) in a certain background. The result bilinear in the external gaugino 
field A is 



(-eV2eA*i;X) x - x = 1 ■ (-*) • (-2e 2 ) 
2^e 2 



d^p X e i>X e 

(2TT) i P 2 (p+kf 



(4tt) ; 



-X e ftX e . 



(29) 



where A e is understood as Fourier transformation (with momentum suppressed) of the external field A e (x). We have 
momentum k flow into the vertex — e\* r 2eA*ipX and for simplicity, we set the momentum exchange through the external 
field A e to be and that through A e to be k. Similarly we have scalar contribution, 



(-eV2eA*ipX) A A* = i-i - (~2e 2 ) 

U 



d 4 p Tr[tf +fljf\ 



(2tt) 4 p 2 {p + k) 2 



A e (A e y 



(47lf 



e 2 k A A e (A e Y 



(30) 
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Again, the momentum exchange through (A e )* is set to be for simplicity. Note that the expectation value propor- 
tional to ■0 e 'i/5 e is the same as Eq. ([2T))) with the replacement of A e — > i\) e . Combine these results, the total expectation 
value is 

9?7> 2 Oj'p 2 4j 

(-eV2eA^X) = j^X'fr' + j^y 2 W + j^^A^f 



2ie 



(X e fX e - \ip e f(ip e + \k 2 A e {A e )*) + ^ e e + k 2 A e {A e )*) (:-!] ) 



(4tt) 2 

In the convention we are using, an operator with a momentum inflow k is given by 



-ik-x 



x O{x). 



(2^)4 

A current with the same charge +1 assigned to ip and A (i.e., Konishi current) is 

iA d A* + ipcr^Tp -> -ik 2 AA* - iifjfiip. 

So it is clear that Eq.([3"Tj) is a linear combination of i? M and the Konishi current Rj[. 

Eq. (j3"Tj) is also consistent with Eq. (fT2")) , at least loosely. So the overall contribution to Jft[{<5^ (<r M J^))] is 

"\p 2p 2 
(e^=A>A) -> (32) 

On the other hand, Eq. (fl"2|) predicts a correction to the U(l) current R'^, 



Pp^-TrW* WlfD 2 {*e v <$>) -+ 99\- ^Jf 2 + ^GMO + 



The first term is part of the trace anomaly Compare this with the superconformal algebra Eq. (|25p . we know 



that the correction to R' is 



which agrees with Eq. (|32|) . So in some sense, we correctly calculate the — | r )fD 2 {t&e v <&) term, though only in a 
way that is not manifestly supersymmetric. 

However, our calculation is not without flaw. We use — ier''(x/j + o il a v Xv) with \n from W a a, the gauge field part 
of the supercurrent to get contribution to the supersymmetry current from the gauge sector while the contribution 
from the matter sector is obtained from a direct generalization of Eq. (|A8I) to 4 + e dimension. Without this double 
standard, the coefficient in front of d^fycr^ip) in Eq. (|3"2"|) will be different and does not match that predicted by 
Eq. (1121) . The physical implication remains the same; namely that R' receives correction proportional to the Konishi 
current and is no longer the i?-current that satisfies the Adler-Bardeen theorem. 

It is not clear whether the operator defined by Eq.tfTUf generates a new i?-symmetry at the fixed point. The 
terminology we use may be a little confusing. By saying "new i?-symmetry" , we refer to the U(l) symmetry that 
forms the superconformal algebra together with supersymmetry and the scaling, instead of the one that transforms 
the gaugino and matter fields according to the charge ratios of 1 : | : — |. The latter is referred to as i?-symmetry. 
Anyway, it is not out of question that the operator of Eq. (fTU)) can be the right current to generate the new i?-symmetry. 
Note that at the fixed point, the trace 0^ scales various fields according to their "quantum dimensions" instead of 
their canonical dimensions though the operator form of this dilatation current is defined according to the canonical 
dimensions. 

In any explicit calculation, it is hard to see how the charges associated with the operator defined by Eq. (flO|) can 
receive quantum corrections. So this operator, after some renormalization, is likely to be the i? M that generates the 
anomalous chiral U(l) symmetry as it is the case in QED. Anyway, in our opinion, the point is that there should be 
an anomalous current that transforms the fields according to the charge 1 : | : — | and it satisfies the Adler-Bardeen 
theorem. Moreover, the latter is definitely not in the same multiplet as the stress tensor. The former may or may not 
be generated by the bare operator defined by Eq. (fTU)) but this could depend on the calculation scheme and is hardly 
physically relevant. 
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A. A Manifestly Supersymmetric Derivation 

The calculation above can be done using the dimensional reduction and the background field method. In [l3j], the 
anomalous dimension term is argued to be zero because of the assumption of on-shell external fields. The assumption 
of on-shell external fields is in general not justified and in this particular case, leads to the missing of a term that 
has physical meaning. In this subsection we recover the anomalous term so that the super-anomaly equation of the 
current is exactly of the form of Eq. (fT3")) . 

In this subsection, we are going to use the convention in [131 ] ■ So we need to determine the corresponding form (in 
this new convention) of Eq.(fT3]). According to Eq.(g5]), and Eq.(C39) in [T3 |. 

V 2 [f • f] = -eW a W a . 

we have the Konishi anomaly, 

V 2 (l>$) = -T{R)[W\ (33) 

Let [J^ a a] be the supcrcurrcnt, rcnormalized so that its expectation value is finite and is exactly equal to what 
one would get by putting into J a a the external fields alone 7 . Note that [Jaa] is the supercurrent that contains the 
energy-momentum tensor in its 69 component . Then the trace anomaly is the 6 component of its super-trace, which 
is given by Eq.(3.15) in |13| . 

V«[J a6l ] = \^V a [W 2 ]. (34) 
9 

The matter contribution to the /3-function, up to two-loop, is (from Eq.fTTJ}) 

% = -3T(G) + T(i2)(l- 7 ) + ... ) 

r 

where 7 is the anomalous dimension (defined from the anomalous scaling of the renormalizcd operators) The scaling 
dimension of the renormalized field $ is 1 +7/2. In order for the missing (missing on the rhs of the anomaly equation 
Eq. ({3"4"|) ) term, V 2 [$$] to give the correct contribution to the /3 function following Eq. (|3"3")) (see [H]]), we require 
Eq. (f34]) to be modified as, 



V* [Jac] = tt m V Q [W 2 ] + 7 V Q V 2 [**] . (35) 



3 W 

Note that the numerical factors in this form are slightly different from those in Eq. (|13p . Now both terms on the rhs 
follow from the vev of 

9 

As explained in section III, the anomaly is determined by W a K a a and from 

V"K a6l = \^ a W 2 + ^V 2 V Q $e v $ + 4ie$e v W a $, 
9 2 

we get — p-VcVF 2 . Let us now show the expectation value of the latter can indeed give the correct super- anomaly 
equation (f3"5"j) . The contribution proportional to V 2 i> e( I? e in the vev of W 2 can be obtained in a similar way as that 
of W a a- Now instead of 

w a6l -> (v 2 v [i y)(v 2 v a y) + ..., 

we have the expansion, 

w 2 -> I(v 2 v Q y)(v 2 v Q y) + . . . , 



7 Here, certain non-local contributions are ignored in the previous literature. 
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as the vertex, where V is the quantum fluctuation of the gauge field. The relevant diagrams are 1(d) (from two 
vertices of <&V$) and 1(e) (one vertex of ^&V 2 $) in [13]. It is not hard to see only the latter gives nonvanishing 
contribution. We have, following the Feynman rule, 



that the vev of W 2 is 



(VV) = 2g 2 d-\ 

(w 2 )^ = i(e) = 2 5 2 $ e n- 1 V Q v^v 2 v a n- 1 $ e 

= 45 2 V 2 $ e n" 1 V 2 V 2 D _1 <I> e 



= -C{R)g 2 V 2 <$> e $ e (36) 
e 

In other words, the last term on the rhs of Eq. (|3l)l) is supposed to be 

(-4v Q W 2 }$$ = -AC(R)g 2 V a V 2 &<S> e . (37) 
9 

The anomalous dimension 7 is given in [l4| as, 

1 = -C(R)-^-AC(R)g 2 , 

where we recall that a = g 2 / (47r). In the last step, we use the convention (47r) 2 = 1 in [l3[. One can see that Eq. ([57| 
exactly agrees with Eq. (|3"5|) . 

B. Charges at the Infrared Fixed Point 

Previously, we have shown that the current R' is a linear combination of the i?-current R^ and R^. In this 
subsection, we apply our result to study the current R' at the infrared fixed point of an SU(N) SYM that has Nf 

matter fields Qf in the fundamental representation and Nf matter fields Q f in the anti- fundamental representation 
8 . The current R' is shown to be the anomaly-free current, whose charge ratios for A, A, ip is, 

l:»l^L:->L. (38) 

We then argue for the advantage of our method compared to the argument [26| based on the approach in [l4| . 
At the infrared fixed point, we have the current R' as, 

R'„ = R» + pRf. 



This follows from the coefficient 3T(G) — 2/(1 — lf)T{Rf) in Eq. (fTB|) (and Eq. ([H|l ). For later convenience, let us 
rewrite Eq. ()13|) in the form, 

qqq n- _ _ i n 2 f _ 1 f 3T(G)-X] / (l-7/)T(_R / ) ] , 2 , 
V Jaa — 3"^ \ TE^y 1-T(G)q/2tt J \- X1 VY J 

-\Y,flfD 2 [{&e v <S> f )]}. (39) 

From the lhs of Eq. (l3T)l) . we have an operator Id^R'^ 9 . Taking its expectation value, the lowest order term is 2d^R^ 
(i.e., Eq. (fTU|) with all fields replaced by their external counterparts). In the context of Slavnov- Taylor identity in the 
background field method [2j, this term corresponds to the contact term and tells us the -R-charges. Moreover, although 
the Adler-Bardeen theorem does not hold for this current R'^, the lowest component of Eq. ([39|) still gives the chiral 
anomaly equation up to one-loop. So a connection between the factor in front of d^R^ and the coefficient of [TrIF 2 ] 



° For a review on supersymmetric QCD and especially the properties at the infrared non-trivial fixed points, see e.g. [29 
9 In fact, we don't really need to know the factor in front of d^R' . 
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can still be established. From Eq.© and the ratios of 3T(G)/jT(Rf) in Eq. (l39l) . we can infer that there is another 
renormalized operator ^[d^Rf] coming out of the term — | jfD 2 [(<&fe v $f)]. The reason is that assigns 

charge +1 to A while R^ assigns charge +1 to ipf and a combination of R^ + ^7-R^ gives the correct coefficient in the 
chiral anomaly equation (coefficient of [TrPF 2 ]). Note that according to Eq. ([TTT) the vev of the operator D 2 [§f e v f) 
is going to give a term proportional to Tr[PF 2 ] and a term proportional to D 2 {(<&fe v Qf)}. The 99 component of these 
two superhelds have f 2 and d^R^ respectively with the appropriate coefficient determined by Eq.©. 

Anyway, at the fixed point, there is no other contribution (from the anomaly) of the form of d aa (\ a \a) 10 and we 
only have d^R^ + ^ r yd tl {R^ 1 ) e (up to a factor) in the vev of D a J a a- As a result, the i?'-charge at the fixed point for 
<f> is |(1 + ^7). At a different scales, the i?'-charge is different. Naively this operator R'^ behaves as a current that 
has different charges when it acts on states of different energy scales. 

With anomalous dimension at the fixed point being, 

37V 

which is necessary for the NSVZ /3-function to vanish, the charge of $ is, 

2, 1 . N f -N 
3 (1+ 2^ = ^VT- 

This result agrees with Eg. (13"5)) . 

The -R'-charges at the fixed point can be obtained in a very different way (2r3 |. In this case, a conserved current is 
defined for every different 7 (see Eq. (2.114) in [HI), 

f q 3T(G)-Z f (l-7f)T(Rf) 

At the infrared fixed point, the second term vanishes and this current is just the supercurrent J aa . However, the 
i?'-charges are obtained from the form of J aa at the UV fixed point, where = 0. It is not clear why this works 
because J a a with different 7/ are different operators. In other words, i?'-symmetries for different J aa arc different. 
Why the charges of one symmetry is determined by that of the other needs to be explained. In our approach, the 
values of the charges for R'^ come out naturally. 

V. PROBLEMS AND CORRESPONDING SOLUTIONS IN THE APPROACH USING TWO 

SUPERCURRENTS 

In this section, we justify the use of the equation of motion in [l2j and [l3| . An equation of motion, when inserted 
into the n-point functions serves like a functional derivative. For example, the expectation value of the equation of 
motion of a field <j>, denoted as $ satisfies 

<^(*W = <^<*>, 

where X denotes other operators. In the background field method, the expectation value of the equation of motion 
of <j) gives something like 

l_ *^T' (40) 

where T is the effective action. In the standard non-supersymmetric calculation of the chiral anomaly using dimensional 
regularization, the chiral current is no longer conserved. In other word, d^jn 7^ by the equation of motion. Instead, 
we have (ip being a Dirac spinor), 

= ~h 5 P^ + H.c. + ^{^, 7 5 }V 



As explained above, the anomaly term [TrW 2 ] has a contribution to the gaugino U(l) current. 
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The first term on the right hand side and its Hermitian conjugate are proportional to the equation of motion and 
both vanish on-shell. The insertion of these terms in a Green's function only gives contact terms (or in the context 
of expectation value, only trivial terms like Eq. ([4"0)l ) but not any anomalous contributions. 

However, in [l2| and (l3j expectation values of the equations of motion apparently do have contributions from the 
anomaly terms. This can be seen as follows. From Eq. (l23l) . we have 



(V<W Q(i ) ~ T{G)W 2 eyit ± 0. (41) 

The right hand side is the anomalous contribution (given by the external fields). Note that assuming no matter fields, 
the equation of motion of the gauge field implies, 

V & W a<i = 0. (42) 

In [l2j and [lH, the lowest order contribution to the expectation value vanishes because of the on-shell assumption 
and a trivial expectation value (no anomaly) means a vanishing expectation value. Therefore, one expects 

(V d ^ Qd ) = 0, (43) 

which is in contradiction with Eq. (|4ip . Moreover, when the equation of motion is used assuming that the expectation 
value of the equation of motion is trivial, which is necessary in the calculation of anomaly, the nontrivial expectation 
value leads to an inconsistency. For example, Eq. ([42]) is used on an operator 

V d J a& = V d (w a& - ^K a& 

to get — |V Q W 2 (coming from the second term — ^W a K a a), which is then taken to be the expectation value of V" J a &. 
In other words, Eq. (|43| is assumed, which is apparently inconsistent with Eq. (j4Tj) . 

In fact, in the scheme of dimensional reduction, the U(1)r current is actually conserved (by the equation of motion) 
and one may expect trivial expectation values for d^R^ and the supertrace of the supercurrent that has i? M as its 
lowest component. This is possible only when we consider the non-local contributions to the expectation values. 

In [281 ] . non-local contributions to the expectation values are considered but do not give any divergent contributions. 
The point is that it is necessary to include the non-local contributions in this calculation scheme (background field 
method and dimensional reduction). For simplicity, the chiral anomaly in (non-supersymmetric) QED is considered, 
as it was the same example that was discussed in [l^ . The non-local contribution to the expectation value of the 
chiral current ipatpa can be evaluated (in this calculation the superscript e for the external field is dropped), 

= -V* /// ^/^(^r^-^^'^VVi 

2dy / dz d i qd i k(-iyq-izk) l3 afi3 1 (q)F a (k-q) 



o Jo 



x / d 4 l- * ~- lkx 



1 

e 



' (I 2 - A) 3 

1 f\dy (""dz ((^kt^M^LJL^ 



9 / I / / -y^b 2 

Zi Jq Jo J J 4> I\j -on, 

= _£ J I d 4 qd 4 k lk ^^M)P"^ ~ g) e -«fcx (44) 

where A = y 2 q 2 — yq 2 + z 2 k 2 — zk 2 + 2yzq ■ k. In the sixth line, we remove the q dependence. Anyway, in the position 
space, the non-local contribution reads, 

Note that the Ad term, which we did not consider, has contribution too. But it appears that this contribution is 
proportional to k a aF 2 , so after a procedure to make (ipaifi a ) real, this contribution drops. 
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The non-local contribution is opposite to the contribution from the e-dimension operators, which is given in [12j. 
So an operator equation like, 

is valid. A renormalized operator [j^] defined in [l2j has a nonvanishing expectation value that gives the correct chiral 
anomaly. 

This result can be generalized to supersymmetric theories. In principle, one can compute the non-local contributions 
to the supercurrent and show that its supertrace vanishes. Such contributions to some of the operators have actually 
been worked out in the literature. For example, the expectation value of the operator $e v $ has a (one- loop) non-local 
contribution (see e.g. Eq.(6.7.10) in 27]). The contribution from the e-dimension operator can be found in Eq.(A28) 
in [H» 

(f$r)(i) = -2C(i?)-U e $ e - ^T(R)-t e • f e , (45) 

These two contributions are opposite to each other. So we have 

£> 2 ($e y $) = 0. 
On the other hand, a renormalized operator [$e 1/ <I > ], 

[i>e v $] = $e y $ - -T(R)T ■ f. 

e 

can be defined to provide the correct Konishi anomaly Eq. (fTTj) . Here we use a bare field $ and therefore the first 
term on the right of Eq. (|45"T) in Ujj is removed. 

The calculations of the non-local contribution that have be given in the literature were not without flaw. The 
non-local contribution to (<&e v <&) is given by the same graphs that contributes to the effective action (Eq. (8) in [28[ , 
more explicitly Ii,^,!^). In [23], it appears that only part of the contribution (1$) is considered. Moreover, I2 is 
infrared divergent. It is not clear whether one can just drop this term. 

However, non-local contributions to the expectation values of operators are infrared finite because the 4-momentum 
injected into the operators becomes an infrared cutoff. As shown by explicit calculation, the expectation value of 
operator $e v $ has terms containing W a W a (W a being the background field). These terms are finite and remain 
non-local after a differentiation by D 2 . We would like to see whether they have an impact on the anomaly. In fl2l |. 
they are discarded with the assumption of on-shell external fields. 

Let us take a look at the 99 component of VW^V^W^. Note that all the fields are functions of = x^ L + i9a^9. 
The expansion around x gives the 9 2 9 2 component. These two components (99, 9 2 9 2 ) are all that are relevant in the 
chiral anomaly equation. First of all, let us take a look at the bosonic contribution, which is of the form of df 2 , 
from V a W tt V^H^. To get a nonvanishing 99 component, we need one 9 and one 9. The lowest component of V Q W Q 
vanishes identically. Therefore both 9 and 9 can not come from one of the V a W a . The 9 (or 9) component of any 
scalar superfield like V Q W Q is certainly fermionic. So there is no bosonic contribution to the 99 component (and the 
9 2 9 2 component) of V a W a V^Wp. A similar argument can be made for (V a V^W]g)W a . 

In summary, we found that the terms containing V a W Q do not give contributions to the chiral anomaly. Nor do 
they give any contributions proportional to F 2 , which would appear in the trace anomaly. Note that the derivative 
expansion does not apply to those V a W a terms because of the infrared divergence in the expansion coefficients. 
As a result those W a W a terms don't give local contributions even after being acted on by D 2 . These terms have 
contributions quadratic in the spinor fields. It is not clear what their physical meanings are. Naively there can be such 
contributions to the one-loop expectation value of the chiral current because of the Yukawa coupling. Even though 
these contributions are just contact terms, the equation of motion is still spoiled. So in order to use the dimensional 
reduction, it appears that we have to assume that the external fields arc on-shell, satisfying the classical equations of 
motion. 

Note that other relevant contributions to the expectation value of $e y $ are proportional to V^W^'V r a Wp). Let 
us also consider the 99 component of V^W^V ( a W^, which is 

V^W ?) V {a W n -► f ap i9«d a& (9y i0 ) cx i9 a 9 & d a& f 2 . 

So we have 

V^W^V^W^+h.c^ 9 a 9°d a& (if 2 -if 2 ). 
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With the assumption of D a W a = 0, V^W^ V (a W p) can be expressed as D 2 W 2 and we can get to the usual form of 
the Konishi anomaly Eq. flTTT) . 

Despite these technical difficulties, it is quite tempting to expect similar non-local contributions to the expectation 
values of W a a , P a a and Q a a ■ These non-local contributions cancel the contributions proportional to K a a and make 
the use of the equations of motion justified. 

VI. RESCALING ANOMALY AS AN INFRARED EFFECT 

The /3-function obtained by the method introduced in [lij agrees with the one obtained by a different method in 
(ltSj . Moreover, the role played by the field strength renormalization is almost the same. Take SQED as an example. 
As pointed out in [l6j the rescaling of the field strength will change the coupling constant and the new theory (at 
cutoff A) is canonically normalized, i.e., has Z = 1. For such a theory with Z — 1 both calculations predict that 
there is no quantum correction for process at the scale A. The effective Lagrangian after the rescaling is exactly the 
effective action for the external fields at the scale A. 

The calculation (of the Jacobian) in [l6[ involves UV regularization and the infrared effects seem to be irrelevant. 
However, we will show that if the momentum modes below an arbitrary scale A are ignored, there will be no con- 
tribution to the Jacobian except for some non-renormalizable terms. The idea is to separate the contribution to the 
Jacobian from the modes above A. To do this, one can consider two Jacobians under the rescaling of field strengths 
at cutoffs A and A' (assuming A < A') respectively. The Jacobian for the scaling of a chiral superfield $/ (in a 
representation of i?/) can be computed (Eq.(A.20) in [lH), 

logJ(A,e"( A >) =-i J d 2 6?W-log(e^)W 2 + 0(±). 

Note that J(A, e Q ( A )) can be understood as the Jacobian from the rescaling (by a factor of e Q< - A ^) of the momentum 
modes k < A. We can set the scaling factor to be the same, then J(A, e"( A ') is the contribution to J(A', e"' A ') by 
modes k < A. The difference log J(A',e"( A — log J(A, e Q ' A ■*) starts from the — -j^, which is the difference of 
the coefficient of a certain non-renormalizable term. This is consistent with the result about the contributions from 
the modes between A and A' in the Wilsonian renormalization group flow, in the sense that the effective theory with 
a cutoff A generally has non-renormalizable terms proportional to the negative powers of A. 

The point is that the multi-loop contributions proportional to log ^ do not come from the rescaling of modes 
between A and A'. The choice of A is arbitrary. In other words, the multi-loop contributions to the /3-function do not 
come from any modes k > A. Therefore, the result in [l6| must also come from the infrared modes as in [l4| 

VII. CONCLUSION 

In this paper, we have tried to elucidate and settle three problems that are related to the anomaly puzzle in J\f = 1 
SYM. First, we study the properties of the current operator R'^ that is in the same super-multiplet as the stress tensor. 
We show explicitly that R 1 ^ is not the same as the (anomalous) current i? M which transforms the fields according to the 
charge ratios 1 : | : — |. Only the anomaly of the latter current is of one loop order and satisfies the Adler-Bardeen 
theorem, while the anomaly of R'^ is proportional to the /3 function. By explicit calculation, we show that R'^ is a 
mixing of the i?-current, and the Konishi current. Moreover, we show that the — | J^f 7/£> 2 ($e^$) term that 
appears in the anomaly equation in [14j gives the same mixed current R'^ and therefore supports the existence of 
two different "supercurrents," even though only one supercurrent was proposed in [l4j |. We then use supersymmetric 
QCD at the infrared fixed point, as an example, to show how the difference between R'^ and i? M can naturally be 
explained in terms of two-supercurrents. 

Secondly, we show that non-local terms must be included for consistency when using the equations of motion in [ll[ 
and [HI . This is necessary because the equations of motion are used there with the assumption that their expectation 
values trivially vanish, while they actually vanish only when the non-local contributions to the expectation values are 
included. 

Finally, we compared the two different calculations of the NSVZ /3 function in [l3| and [l6j]. The second method, 
which is based on the Jacobian arising from field strength rescaling, seems independent of the infrared behavior of the 
theory, while the first method seems to depend only on the infrared behavior. We resolve this apparent contradiction 
by showing that the infrared modes are also crucial in getting the multi-loop corrections to the /3 function in the 
second method. The reason, as we show, is that the contributions from modes above any arbitrary nonzero scale, 
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A, to the rescaling Jacobian are proportional to non-renormalizable terms and therefore do not contribute to the /3 
function. 
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Appendix A: Supercurrent 

Here we will review some properties of the supercurrent. We mostly use the convention in !19j including the choice of 
matrices and superderivatives. The only two different choices are the form of vector superfield and the integration 
of Grassmann variables. First of all, let us put down the transformation rules for components of a chiral field $ that 
appears in the Wess-Zumino model with the Lagrangian, 

C = \J d 2 6d 2 6<S>$> + i J d 2 9g<$> 3 . (Al) 

The chiral superfield can be written in the component form, 

$ = A + V2(9?A + 6 2 F, (A2) 

where the factor of y/2 is to make ip canonically normalized. The supersymmetry transformation (parameterized by 
of the components is given by (g — for simplicity), 

SA = V2^ip, 5 A* = (A3) 



Sip = iV2a"^A + V2£F, Sip = -i^a^d^A* + V2£F* 



(A4) 



SF = iV2^d„ip, SF* = iV2^df,ip 
With this transformation, we can work out the current from Lagrangian, 

C = id n ipa n ijj + A*aA + F*F 

= ^m^V " y>&*W - d„A*d^A + F*F, 

which is just Eq. (|Al[) in component fields. 

For simplicity, we set F — F* = 0. Consider only the transformation generated by £, we have 



<5f£ = a/2 
= V2 
= V2 

So the current is given by, 



d^A*d»m + l -d^-iio»d v A*)d^ - l -(-i^d v A*)d^ip 
d^A*d»m + ^{d^)a v d v A*a^ + ^d^A*^ - ^d^cr" d v A*a^) 



J„ = sf2d v A*a v aa<<P- 



(A5) 



(A6) 



(A7) 



For a general theory, we need to keep the auxiliary fields when we do the variation. For a free theory, we can see that 
Eq. (IA6l) implies no F in J. The variation of the kinetic term ipa^d^ip contains something proportional to d^F). 
The part with (d^F)^ is combined with another term to form a total derivative d^K^ but this does not appear in 
the conserved current since it has to be the same as the part (from the variation of ipa^d^ip) proportional to (d^)F. 
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Despite that, we can use equation of motion in the interacting theory to get F in J M , the current following Noether 
method does not contain F . All the interacting terms containing F do not have derivatives and therefore cancel 
identically. The extra terms (proportional to coupling constant) in J M involves variation of F though. 

Note that the charge generated by this current Q = J J° does generate the correct supersymmetry transformation 
on component fields. This is obvious for ip and A. One subtlety is the transformation on A. It gives 

5(A) = -jv / 2£[<9 i AVffo</', A] = V^tao^diip = -V2£aoa°d t 2p = 



where we have used the equation of motion. 
However, we actually use the J p given by, 

J a = V2 



d v A*a v a^ + ^ v d v {A*i>) 



(A8) 



This is the "improved" supersymmetry current of which gives the same charge and is also conserved. Note that 
the second term does not have any contribution to the supersymmetry charge Q because we have a 00 = 0. Although 
<t 7^ 0, the spatial derivatives do not contribute to the charge since they only give boundary terms when integrated 
over a spatial slice (volume). 
The supercurrent is given by, 



Ju 



Operator Eq. (|A9p can certainly be expressed in the form of Eq.Q. Explicitly, we have 

C ■ 



2i A^t A* 2 , 7 

—A d a& A + -lAaVd 



(A9) 



(A10) 



This indicates a correct 2 : — 1 charge ratios for bosonic and fermionic fields. The lowest component is related to 
R-current, 



By the way, with the introduction of both a mass term m<E> 2 and a cubic term g$> 3 , the U(l) transformation is 
no longer a symmetry. A charge given by the above does not give another supercharge Q when we take their 
commutator [R, Q]. 

The 8 component of this current is not the supersymmetry current J M . Instead, they are related by, 



(All) 



Now let us show this. Both second terms in Eq. (|A8p and Eq. dAlip give currents that are in the same equivalence 
class (second term on the right hand side) as Eq. (|A7[) . 



3 
2 
3 

2a/2 r 



i\/2> 0„ A* - -^=d u A*a u a„tp + V2e afs F^ 

iV2^ dp A* - -^=d u A*a v a^ - ^a^F 

.1 



id^A* - ii-d^a^ip + d u A*a v a^) 



(A12) 



Note that we use D = dg, D = —da — i28d. Now we have 



2\/2 .[ 
— - — i 



(id^A* + ia„a u d v ^A*) - z(-^AV M a> + d v A*a u a^) 



V2 



d v A*a v a^ + ^a^d v (ipA*) 



(A13) 
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which is exactly Eq. (|A8[) . In the third line, we use a p a p = —4 and overdo- 1 ' = 2er M . 

The commutators of Q with the component fields of $ give the supersymmetry transformation of the fields when 
the equation of motion holds. 

The M Qt j for this supercurrent in fact vanishes. The 9 2 component of $ vanishes (F = 0) and therefore the first 
term in J p does not give any contribution. On the other hand, D$ has neither 9 2 nor 9 2 component. So to get a 9 2 , 
we need a 9 from each factor (£>$ and -D$), _D$ only has a nonvanishing 9 component. The 9 component is just F 
and therefore vanishes. Similarly, the 9 component of D<& vanishes and we have a vanishing M^. 

The bosonic part of the 99 component of the supercurrent can be worked out following from Eq. (|A9[) . 



2d a$ AS^d p& A* - e^(d g6 A d a6t A* -A d a& 8 gB A*) 



fsfs J 



(A14) 



We use v^u to denote the 99 component and t pv = hvi^u) a s one half of the symmetric part of the 99 component. 
With a tpv defined as 



- 2d {tl Ad v) A* + \v^d p Ad p A* + ^d v AA* + ^Ad„d„A* 



(A15) 



we can obtain the 66 component explicitly, 



% % 



daaA dppA — dppA d a aA — 2d a pAdfiaA + d a adppA A + Ad a adppA 



(A16) 



Note that this agrees with Eq. (|A14|) . 

Let take a look at the relationship between t pv and the stress tensor. The stress tensor can be obtained as, 



Tpy = ~d^Ad v A* - d v Ad p A* + r lpu 8 a Ad a A* - -(^8^ + ^d p ^) 

% — — % — % — 

--(ipa^ip + ipa u dpip) + ^^{-ip^dpip + -ipa p d p ip). 



(A17) 



Of course, we usually use 'duu, which is defined as the improved stress tensor of T, 



////. 



<V = T pu + ^(d^- V ^D)AA* 

= ~d^Ad v A* - 2 -d v Ad p A* + ^ ntlu d p Ad p A* 

+±d tt d v AA* + ^Ad p d u A* - ^(ADA* + DA A*). 



(A18) 



The extra term to T pu does not provide any contribution to the charges P p and certainly does not affect the conser- 
vation. Note that we have t p p = \(ADA* + DA A*), and therefore 



V - Wp P = ~^9 {fl Ad v) A* + -r ]fW d p Ad p A* + -d p d„AA* + -Ad p d„A* 



--^{ADA* +DAA* 



= 0„ 



(A19) 



Of course, following Eq. (|A19[ ). we have proportional to ADA* + DAA*. Similar conclusion should hold for 
improved supersymmetry current J p . 

Note that this difference between t M1/ and follows from the difference between x M and J p Eq. dAlip . This 
statement is actually only true when is symmetric. With Eq. (IA19p the variation of J M (under a supersymmetry 
transformation generated by £) is proportional to while that of \ p is proportional to 

Let us now include the fermionic fields in our considerations. The 99 component, v pv (9a v 9\ of the supercurrent is 
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(following from Eq. (|A9j> ). 



V UU W(7 



- i2 ■ 2{9a v ) a d v {9^){(j^) a + i2 ■ 2{^a ll ) & {9a v ) a d v {9i)) 

2 r 

»3 



— % — % — % — % — 

iip(j v dpi/) — -d^ipa^ + -ipa^dvi/j - -d p \j}a v a p <7 tl \jj - -d p ^a„(T p a^ 



- 2 ( _ j i 

1 - i - 

-^vp P ad P ^O a i) - ~d p 1p(r) vfl (T P - 5 P (T V - + H . C . 

We would like to show that the anti-symmetric part of v^v satisfies the following expression on-shell, 



1 



v [nu\ — ^tJ.uabd a R b 

It is not hard to rewrite he^ a bd a R b in the form given by, 



(A20) 



(A21) 



Following, 



we have 



e f ivabd a R — -e a /30 7 0-R 7< i) + H. c . 



<TL*;^a b d a R b = -^(a^eUe. $ + (ean^ a0 ]e, uab d a R b 



Now we have 



= «e Q /39 7 ( fi (?A 7 V5 / 3 ) ) - id 1 (aR^pfafl 



(A22) 



(A23) 



(A24) 



This is the part antisymmetric in a, and symmetric in a, j3 and corresponds to the self-imaginary dual component 
of the t^u] ■ The Hermitian conjugate gives another term, 

-2* £ A/j 5 (" 7i? /3)7 



With the use of Eq.(TA24| and Eq.(|J22j|, we get Eq.flHI]). 

Note that in the derivation above, we did use equation of motion. In fact, Eq. (|A2ip is true on-shell. It follows from 



which can be derived using the explicit form Eq. (|A9[) and equation of motion (and its conjugate) 

L> 2 $ = 0, 

and, 

[Da,D 2 ] — AiD a a l ^ t6l d ll . 



(A25) 
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In component (more explicitly 6 component), Eq. (IA25P implies 

i{a u 6) a d u R a& + {a v 6) a v a&v = iduR p eed v (j' 1 + v^Sea'a" 

= d a R b e ab ^8(e<j^) + 2v [ ^ ] 9(ea^) 



0. 



(A26) 



This is equivalent to 



which is exactly Eq. (|A2ip . 

Let us now consider SQED, whose supercurrent is given by 



We have R a a = — -§iX a Xa- The 96 component is decomposed in a way as Eq.@ (with only fermionic fields considered), 
26 f3 9^(iX a d g g\a ~ idgg\ a \a) = 26 p 9^\(iX( a d g) gXa-id g( 6X a X 6( -\) 



1 



where r M „ is given by, 



4 " " e 2 



29"^ + -eJ^d^R^+H.c), 



(iX a d g gXa + iXpd a pXa — idgpX a Xa — idfiaXaXp) 



-{^apd^pX^Xa+K.c) 

1 - i 1 i 

-—ziiXdvO-^X - -Xa v a p (j^d p X + H.c) - -^{-d p Xa p u u a IJ ,X + H. c) 

1 % — % — 1 % — % — 

~^1^ dv<T i Ji ^~ ■^Xa l/ (j p a^d p X) - -^{-Xa l _ l a p a u d p X+ -Xa^X) +H.c 
-^C-Xd^^X + YL.c) - ^ €pvba d b {-^Xcr a X). 



(A27) 



(A28) 



Note that we have a — je /1 ^ a bd a R b in t^ u . So t mi , is not symmetric. Following a similar derivation as in the Wess- 
Zumino model, we can show that the tui V ~\ — and certainly, the symmetric part produces the stress tensor. This 
symmetric part t^ v = is in fact related to the stress tensor in the same way as Eq. (|Al9j) . 
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